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abstract 

We present and discuss the photometric and spectroscopic evolution of the peculiar SN 1998bw, asso- 
ciated with GRB 980425, through an analysis of optical and near IR data collected at ESO-La Silla. 

The spectroscopic data, spanning the period from day —9 to day +376 (relative to B maximum), have 
shown that this SN was unprecedented, although somewhat similar to SN 1997ef. Maximum expansion 
velocities as high as 3 x 10^ km s~^ to some extent mask its resemblance to other Type Ic SNe. At 
intermediate phases, between photospheric and fully nebular, the expansion velocities 10^ km s~^) 
remained exceptionally high compared to those of other recorded core-collapse SNe at a similar phase. 
The mild linear polarization detected at early epochs suggests the presence of asymmetry in the emitting 
material. The degree of asymmetry, however, cannot be decoded from these measurements alone. 

The He I 1.083 fj,m and 2.058 fj,m lines are identified and He is suggested to lie in an outer region 
of the envelope. The temporal behavior of the fluxes and profiles of emission lines of Mg I] 4571 A, 
[O I] 6300,6364 A and a feature ascribed to Fe are traced to stimulate future modeling work. 

The uniqueness of SN 1998bw became less obvious once it entered the fully nebular phase (after one 
year) when it was very similar to other Type Ib/c-llb objects, such as the Type lb SN 1996N and 
the Type lib SN 1993 J, even though SN 1998bw was 1.4 magnitudes brighter than SN 1993 J and 3 
magnitudes brighter than SN 1996N at a comparable phase. 

The late phase optical photometry, which extends up to 403 days after B maximum, shows that the 
SN luminosity declined exponentially but substantially faster than the decay rate of ^®Co. The UVOIR 
bolometric light curve, constructed using all available optical data and the early JHK photometry 
presented in this work, shows a slight flattening starting on about day -1-300. Since no clear evidence 
of ejecta-wind interaction was found in the late time spectroscopy (see also Sollerman et al. 2000), this 
may be due to the contribution of the positrons as most 7-rays escape thermalization at this phase. A 
contribution from the superposed HII region can, however, not be excluded. 

Subject headings: Supernovae: general; supernovae: individual SN 1998bw; gamma-rays: bursts 



1. INTRODUCTION 

SN 1998bw was discovered by Galama et al. (1998a) 
in the BeppoSAX Wide Field Camera error box of 
GRB 980425 (Soffita et al. 1998, Plan et al. 1999) close to 
a spiral arm of the barred galaxy ESO 184-G82, by com- 
paring two frames taken at the ESO New Technology Tele- 
scope (NTT) on Apr 28.4 and May 1.3 UT. Spectroscopic 
and photometric observations, both in the optical and in 



the near IR, started at ESO-La Silla immediately after 
the discovery, and showed that this object was profoundly 
different from all then known SNe (Lidman et al. 1998). 
The peculiar spectrum led to diverse classifications. A few 
days after its detection, the object was classified as a SN lb 
(Sadler et al. 1998) and later as a peculiar SN Ic (Patat & 
Picmonte 1998a, Filippenko 1998) owing to the complete 
absence of H lines, the weakness of the Si II 6355 A line 
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The metamorphosis of SN 1998bw 



Fig. 1. — Upper panel: SN 1998bw in ESO 184-G82. The frame was obtained on June 17, 1999 (403 days after B maximum light) in the R 
band (5 min) with the ESO 3.6m+EFOSC2. Lower panel; HST+STIS image of SN 1998bw environment at about 764 days after B maximum 
light (Fynbo et al. 2000). The SN location is marked by a small circle and the big circle centered on the SN position has a diameter of 1". 



and no clear He I detection in the optical spectra. 

Its peculiar spectroscopic appearance (Galama ct al. 
1998b, Iwamoto et al. 1998), its unusually high radio 
luminosity at early phases (Kulkarni ct al. 1998) and, 
in particular, the probable association with GRB 980425 
through positional and temporal coincidence (Galama et 
al. 1998b, Plan ct al. 1999) placed SN 1998bw at the 
center of discussion concerning the nature of Gamma Ray 
Bursts (Wheeler 2001). 

Independent photometric and spectroscopic data sets 
have been presented by several authors, whose results 
will be discussed and compared with those presented here 
throughout the paper. Here we give only a brief account 
of the main results of these previous works. 

The early light curve of SN 1998bw has shown that 
the object was unusually bright when compared to known 
SNc of Type Ib/c (My ~ -19.2+51og hes, Galama et 
al. 1998b). The broad-band photometric observations 
by McKenzie & Schaefer (1999) taken during the inter- 
mediate phases (47 171 days after maximum brightness) 
revealed that the object settled on an exponential decay 
similar to that observed in other Ic. McKenzie & Schae- 
fer first suggested that even in this case the light curve 
was powered by the radioactive decay of ^^Co with some 
leakage of 7-rays. Finally, the late phase light curve was 
covered up to 500 days from the explosion by the observa- 
tions of SoUerman et al. (2000). Their models achieved a 
fairly good reproduction of the data with the radioactive 
material well mixed in the ejecta and M(^^Ni)~0.5 M0. 

The peculiar spectroscopic behavior of SN 1998bw 
around maximum light has been presented and discussed 
by Iwamoto et al. (1998), who identified the main spectral 
features as O I, Ca II, Si II and Fe II. The estimated expan- 
sion velocities were exceptionally high ('^SO.OOO km s^^) 
and this caused a severe line blending. The evolution 
during the first months was unusually slow compared to 
known Ic, with the nebular spectra still retaining many of 
the features present during the photospheric phase (Patat 
et al. 2000, Stathakis et al. 2000). The late onset of 
the fully nebular phase has been interpreted as an indica- 
tion for a large ejected mass (Stathakis et al. 2000) as it 
was predicted by the early light curve models (see below). 
Dtiring the intermediate phase, the emission lines were def- 
initely broader than in known Type Ib/c supernovae and 
the simultaneous presence of iron-peak and a-elements in- 
dicated unusual relative abundances or physical conditions 
in the SN ejecta (Patat & Piemonte 1998b, Patat et al. 
2000). The late spectroscopy presented by SoUerman et 
al. (2000) showed that the tentative morphological classi- 
fication of SN 1998bw as a Type Ic event was indeed ap- 
propriate. The main features have been identified as [O I], 
Ca II, Mg I and Na I D, the latter possibly contaminated 
by He I 5876 A. 

Photometric and spectroscopic modeling of the early 
phases and the possible connection with GRB 980425 have 
been discussed by Galama et al. (1998b), Iwamoto et al. 
(1998), Hoflich, Wheeler & Wang (1999), Woosley, East- 

^ IRAF is distributed by the National Optical Astronomy Observat( 



man & Schmidt (1999) and Plan et al. (1999). 

The symmetric models of the early spectra and light 
curve by Iwamoto et al. (1998) and Woosley et al. (1999) 
reached a similar conclusion: SN 1998bw was generated by 
an extremely energetic explosion of a ~10 M© C-0 star, 
which ejected 0.5-0.7 Mq of ^^Ni. This large amount of ra- 
dioactive material; which is one order of magnitude larger 
than that estimated for known core-collapse SNe (Patat 
et al. 1994, Schmidt et al. 1994), is required to power the 
early light curve. In these models the explosion energy 
must be large to accelerate the ejected mass (~10 Mq) to 
the observed velocities and to make the light curve peak 
in only 19 days. Actually, in spherical symmetry, the light 
curve around maximum can be reproduced using different 
combinations of explosion energy and envelope mass. This 
degeneracy was resolved by Iwamoto et al. (1998) using 
velocity information and computing synthetic spectra for 
the various candidate models. 

An alternative scenario was proposed by Hoflich et al. 
(1999), who used less energetic asymmetric models and a 
suitable combination of viewing angle and degree of asym- 
metry. The required mass of ''^Ni is in this case 0.2 Mq. 
This low mass of radioactive material, however, seems too 
low to explain the late emission of the supernova (SoUer- 
man et al. 2000). 

In this paper we report on the results of an extensive ob- 
servational campaign carried out at ESO-La Silla, which 
covered the evolution of SN 1998bw from the discovery 
until 417 days after the Gamma Ray Burst detection. The 
paper is organized as follows. 

In Sec. 2 we discuss the observations and reduction tech- 
niques for the optical and the near IR. The evolution of 
SN 1998bw around maximum light is discussed in Sec. 3, 
4 and 5, which deal with low resohition spectroscopy, po- 
larimetry and high resolution spectroscopy respectively. 
The He I detection in the near IR spectra is presented in 
Sec. 6 in which possible alternatives to the He I identifica- 
tion are also investigated. The description of SN 1998bw's 
metamorphosis ends in Sec. 7, where we follow its evolu- 
tion into the nebular phase. The late phase light curves 
are presented and compared with those of other SNe in 
Sec. 8. This section also presents the UVOIR bolometric 
light curve, which is compared to the model of Iwamoto 
et al. (1998). Finally, Sec. 9 summarizes our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Optical Spectroscopy 

Spectroscopic observations were obtained from early 
to late phases using the ESO 1.52m (Boiler & Chivens 
spectrograph), ESQ -Danish 1.54m (DFOSC), ESO-3.6m 
(EF0SC2) and the ESO NTT (EMMI) telescopes. Expo- 
sure times ranged from 10 minutes near maximum light to 
several hours at late phases. The journal of spectroscopic 
observations is shown in Table 1. 

The spectra were reduced using IRAF'^ packages. For 
some of the spectra taken at the ESO-Danish telescope, 

;, which is operated by the Association of Universities for Research in 
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Fig. 2. — The raw spectrum (ADU), position angle (degrees) and linear polarization degree (%) are shown from the EFOSC2 spcctro- 
polarimetry observations at day —7 (upper panel) and day +10 (lower panel). Wavelengths are corrected to the host galaxy rest frame 
(f9a!=2532 km s~^, see Sec. 5). The linear polarization and position angle data have been binned into 200 channel increments to increase 
the signal-to-noise. 

Fig. 3. — Spectroscopic evolution of SN 1998bw from day —9 to day +22. For presentation the spectra have been vertically shifted by 
arbitrary amounts: -7^ (-0.40, log F^), -6=* (-0.60), -2'* (-0.95), -1'* (-1.10), +1^^ (-1.30), +3^ (-1.45), +4'' (-1.60), +6"* (-1.70), 
+11^* (-1.75), +12<^ (-1.90), +13'' (-2.10), +19'' (-2.30), +22<^ (-2.45). Spectra are in the host galaxy rest frame (?;gu;=2532 km s ^, see 
Sec. 5). 



it was not possible to remove the fringing in the red since 
suitable flat fields were not available. This gives rise to the 
high fr(^qiKmcy modulation of the spectra at wavelengths 
longer than 7500 A. 

Particular care was devoted to the extraction of the SN 
spectra to avoid contamination from the host galaxy back- 
ground. Nevertheless, especially in the late phase spec- 
tra, the contribution from an underlying HIT region could 
not be completely eliminated, and thus unresolved nar- 
row lines (Ha, H/3, [O II], [O III]) appear in the reduced 
spectra. We emphasize that these features do not show a 
coherent time evolution, but rather depend on seeing con- 
ditions and slit position. For this reason wc conclude that 
thoy arc not intrinsic to the SN. Finally, there is no evi- 
dence of continuum contamination from the parent galaxy. 

Wavelength calibration was achieved by using arc spec- 
tra from He-Ne-Ar lamps, while the response curves were 
obtained via observations of spectrophotomctric standard 
stars (Oke 1990; Hamuy et al. 1992). The accuracy of 
the absolute flux calibration was finally checked against 
the broad-band photometry and, when necessary, adjusted 
accordingly. 

On May 6, 1998, two high resolution spectra of 
SN 1998bw, covering the region 3750-7650 A, were ob- 
tained at the ESO-NTT using EMMI in the echelle mode, 
with a 7" long and 2'.'0 wide slit (resolution 1 A FWHM 
at 5900 A, i.e. 50 km s^^). Order definition and ex- 
traction, sky subtraction and wavelength rebinning were 
achieved using the IRAF echelle reduction package. The 
orders were merged and the two exposures combined to 
a weighted mean to give the final spectrum which was 
eventually flux calibrated by comparison with a medium 
resolution spectrum of the SN taken the following day (see 
Table 1). A signal-to-noise ratio of about 15 in the region 
of the Nal D doublet was achieved. 

2.2. IR Photometry and Spectroscopy 

Near IR photometry of SN 1998bw was obtained at three 
epochs with SofI at the ESO NTT (see Table 2) through 
the J, H, and K passbands (Bessell & Brett 1988). In 
order to allow for a proper sky subtraction, the target 
was observed using the jittering technique. The reduc- 
tions were performed in a standard way within IRAF and 
the fluxes were estimated via plain aperture photometry, 
since at those early epochs the contribution from the host 
galaxy was negligible. The photometric errors were es- 
timated including the nightly zero point variations, the 
aperture correction and the photon shot noise. 

Near IR spectra of SN 1998bw were taken at the same 
epochs and with the same instrument used for the IR pho- 
tometry. To cover the entire near IR region, two grisms 
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were used: one for the range 0.95-1.64 /zm and one for the 
range 1.53-2.52 /zm. The observations were made with 
the 1''0 slit. The resulting resolving power, measured from 
the comparison xenon spectrum, was A/ A A ~600 for both 
grisms. 

As is standard procedure for IR spectroscopy, the SN 
was observed at two positions along the slit and the tele- 
scope was nodded between these two positions once every 
few minutes. The highly variable night sky is then accu- 
rately removed by combining all spectra in the appropriate 
way. Atmospheric features have been removed by dividing 
the extracted spectra by the spectra of nearby bright stars 
that were observed soon after the SN: HD 135591 (May 
18, Ferryman et al. 1997), Hipparcos 106725 (June 6, 
Ferryman et al. 1997), BS4620 and BS6823 (June 30, Mc- 
Gregory 1995). The stellar spectra often contain weak ab- 
sorption lines from hydrogen (Faschen and Brackett series) 
and helium, which where removed before division. Wave- 
length calibration was achieved using comparison emission 
line spectra of xenon gas lamps and is accurate to 1-2 A. 
The spectra obtained with the two different grisms were 
then combined into a single spectrum. 

Relative flux calibration was achieved by multiplying the 
spectrum by a blackbody curve with a temperature appro- 
priate for the star used to remove the atmospheric features 
(HD 135591, T=30.000 K; Hipparcos 106725, T=5800 
K; BS4620, T=13,000 K, BS6823, T=25,000 K). Abso- 
lute flux calibration was performed by comparison to the 
broadband IR photometry at the same epochs. 

2.3. Polarimetry 

Spcctropolarimetry of SN 1998bw was performed at two 
epochs (1998 May 4 and 1998 May 29) using the ESO 
3.6m equipped with EF0SC2 in polarimetric; mode (Fatat 
1999). A WoUaston prism was used with a focal plane 
mask to isolate the ordinary and extra-ordinary ray spec- 
tra of object and two sky positions; the separation of the 
two spectra is 20". A Half Wave plate was used to obtain 
spectra at four different position angles of 0?0, 22?5, 45?0 
and 67?5. The B300 grism was used, giving a resolution of 
10 A over the wavelength range 3400-7550 A. Table 1 pro- 
vides a journal of the observations. The airmass was large 
at the time of the observations and the spectral region at 
A <4000 A is severely affected by differential atmospheric 
refraction. The HST polarized standard star HD 161056 
(Turnshek et al. 1990) was observed on each night to check 
the polarization and position angle calibration. 

The data were bias subtracted and the ordinary and 
extra-ordinary ray spectra of SN 1998bw were extracted 
with the local sky subtracted. Statistical errors were as- 
signed to each extracted spectrum using the known prop- 
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Fig. 4. — Comparison between SN 1998bw and other SNe at maximum light. Data are from Patat et al. 1995 (1994D), Harkness et al. 
1987 (1984L), ESO-KP data base (19941) and Garnavich ct al. (1997cf, in preparation). All spectra have been correetcd to the parent galaxy 
rest frame adopting the following recession velocities: SN 1998bw, 2532 km s'^ (see Sec. 5); SN 1997ef, 3504 km s'^; SN 19941, 463 km s'^ 
SN 1994D, 448 km s"!; SN 1984L, 1934 km s'^ (De Vaucouleurs et al. 1991). 



Fig. 5.- 
velocity. 



Evolution of the Si II 6355 A region. The empty circles mark the value which has been assumed to represent the photospheric 



Fig. 6.— The Si II 6355 A expansion velocity evolution of SN 1998bw compared with SN 1994D (Patat et al. 1995), SN 19941 (Millard 
et al. 1999) and SN 1997ef (Garnavich et al. in preparation). For comparison the photospheric expansion velocities deduced from the Ha 
absorption minima of SN 1987A are also plotted (Phillips et al. 1988). 



erties of the detector (CCD ESO#40, Patat 1999). The 
data were then combined, following the procedure outlined 
by Tinbergen & Rutten (1992) and the total spectrum, lin- 
ear polarization and position angle were computed using 
dedicated programs running in the ESO MIDAS package 
(see Walsh 1992). The polarization was binned into equal 
wavelength width bins to provide polarization errors per 
bin below 0.1% over the wavelength range 4000-7000 A. 
For observations of HD 161056, the measured polarization 
in the V band was 4.3% at position angle (PA) 76°, in 
satisfactory agreement with the standard value of 4.04% 
at PA=67°, given that Turnshek et al. (1990) suggest it 
may be variable in polarization. 

2.4. Optical Photometry 

Late time optical broad band photometry of SN 1998bw 

in the Johnson-Cousins UBVRI photometric system 
(Bcsscll 1990) was obtained in the phase range 310-403 
days from B maximum. Table 3 shows the journal of 
the photometric observations, which were performed with 
the ESO-Dutch 0.9m telescope, equipped with a TEK 
coated 512 x 512 pixels CCD (scale 0.44'7pixcl), and 
with the ESO-3.6m telescope, equipped with EF0SC2 (Lo- 
ral/Lesser 2048 X 2048 pixels CCD, scale 0.16"/pixel). The 
SN is located on a spiral arm and superimposed on an HII 
region, which is clearly visible in Fig. 1. At late phases the 
background contribution is significant, and the SN magni- 
tude cannot be measured using plain aperture photome- 
try. Therefore, after the standard bias and flat-field correc- 
tions, magnitude measurements were obtained by means 
of the point-spread function fitting technique implemented 
in SNOoPY (Patat 1996), a package specifically designed 
for SN photometry within the IRAF environment. 

Color terms for the two instruments have been com- 
puted using observations of standard fields (Landolt 1992). 
The SN magnitudes have been calibrated by means of rela- 
tive photometry with respect to the local sequence defined 
by Vreeswijk et al. (2001, stars 1 to 10). The results of 
the late photometry here presented are in good agreement 
with those recently published by SoUerman et al. (2000). 

3. SPECTROSCOPIC EVOLUTION AROUND MAXIMUM 

LIGHT 

The spectral behavior of SN 1998bw from day +7 to day 
-|-94 has already been discussed by Stathakis et al. (2000). 
The larger wavelength range and improved phase coverage 
of the data we are presenting here, allow us to analyze the 
early phases in more detail and to push the discussion fur- 
ther into the late evolutionary stages. For a discussion of 

2 The data for SN 19941 are from the ESO Key Project on Supernovae Data Base (Turatto et al. 1990), hereafter ESO-KP. 



an independent late phase data set see SoUerman et al. 
(2000). 

The spectroscopic evolution of SN 1998bw around max- 
imum light is presented in Fig. 3. Wavelengths are cor- 
rected to the parent galaxy rest frame, adopting Vgai = 
2532 km s~^, as measured from the narrow Ha emission 
line (cf. Sec. 5). Phases are from B maximum, whose 
epoch (May 10, 1998 UT) was estimated from the photo- 
metric data presented by Galama et al. (1998b). It oc- 
curred 14.4 days after the detection of GRB 980425. We 
adopt these conventions throughout the paper. 

The general appearance of the spectrum at maximum 
light is quite unique among SNe, even though it is some- 
what reminiscent of SN 1997ef (Filippenko 1997c; Gar- 
navich et al. in preparation) which has been modeled as 
a massive SN Ic event (Iwamoto et al. 2000; Mazzali, 
Iwamoto & Nomoto 2000). This appears clearly in Fig. 4, 
where the spectrum of SN 1998bw is compared to those of 
other Type I SNe. 

Another striking feature is the redward shift with time 
of most of the spectral features, both in absorption and in 
emission. At early phases the apparent broad emission - 
like features do not result from discrete emission lines, but 
occur where absorption line optical depths are low and 
photons red-shifting in the expanding ejecta, as a result 
of absorption and scattering processes, have a higher prob- 
ability of escaping (see for example Mazzali 2000; Pinto & 
Eastman 2000). Thus, while the shift of the absorption 
lines can be understood as the effect of the inward reces- 
sion of the photosphere, that of the emission features re- 
quires a clearer knowledge of which lines are contributing 
to absorption, and how this changes with velocity (i.e. co- 
moving radius). The effect is also shown in Fig. 5, where 
the evolution of the Si II 6355 A region has been plotted. 

In these early phases when the velocity is high, line- 
blending is particularly severe. The modeling presented 
by Iwamoto et al. (1998) suggests that the main features 
are due to lines of Si II, O I, Ca II and Fe II. The time 
evolution of the expansion velocity at the photosphere, as 
deduced from the minimum of the absorption trough of 
the Si II 6355 A line, is shown in Fig. 6. For compari- 
son, the corresponding values measured for the Type la 
SN 1994D, the Type Ic SN 19941^ and the pecufiar Type 
Ic SN 1997ef are also plotted. SN 1998bw shows a sudden 
break in the expansion velocity decline rate around day 
15 and another at day 27. While changes in the slope are 
to be expected as a result of changes in optical depth at 
line center, the reasons for sudden changes are not clear. 
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There is no doubt that the velocity dechne rate changes 
significantly; what is less clear, in view of measurement 
uncertainties, is how suddenly they occurred. 

The velocity deduced from the Si II 6355 A line is about 
30,000 km s~^ at day —7, and decreases to about 8000 
km s~^ at day +22. These values arc exceptionally high, 
for any SN, as can be seen in Fig. 6. SN 19941 shows a 
trend which is very similar to that of SN 1998bw, although 
the velocities are systematically lower, while SN 1997ef 
deviates significantly after day 25 (Fig. 6). Phases for 
SN 1997ef are uncertain because maximum light was not 
recorded. 

4. POLARIMETRY 

In order to investigate any possible asphericity in the ex- 
plosion of SN 1998bw, spectro-polarimetry measurements 
were performed at two epochs (day —7 and day +10, see 
Tab. 1). The polarization averaged over the range 4000- 
7000 A was 0.7% and 0.5% at the two epochs respectively. 
Figure 2 shows the total spectrum (raw, not flux cali- 
brated), the polarization binned into 200 channel bins and 
the position angle of the polarization vector projected on 
the sky at both epochs. The formal errors per bin are 
about 0.04% but the error on the polarization binned over 
the optical range is 0.1% given that the measurements were 
not of high photometric quality and the discrepancy on the 
measurement of the polarized standard. It appears that 
the polarization spectrum at day —7 may have a signifi- 
cant increase to the red which is not in evidence at the 
later epoch. In terms of the integrated polarization over 
the 4000-7000 A range, there may be a marginal change 
in the polarization between the two epochs, but no great 
significance is attached to this result. 

The interstellar polarization in the direction of the host 
galaxy can be estimated by the measurement of stellar 
polarization along a nearby sight line. Kay et al. (1998) 
measured the linear polarization of HD 184100 at 0.75% 
in PA 176.5. Correcting the linear polarization results 
by this value of interstellar polarization results in 0.6% 
at PA 80° and 0.4% at 67° at day -7 and +10 respec- 
tively. The polarization value is not very different from 
the value of 0.53% measured by Kay et al. (1998) at day 
~42 and is most probably intrinsic to the SN, although 
a dusty medium in the SN parent galaxy cannot be ruled 
out. However the position angles differ from the Kay et 
al. value of 49°, although the errors on these earlier deter- 
minations are larger (at least 5°). Despite its association 
with a 7-ray burst, these polarization values are similar 
to those reported for other core-collapse SNc (Wang et al. 
1996). 

The small degree of polarization at optical wavelengths 
can be explained in terms of a moderate departure from 
sphericity (axial ratio <2:1, Hofiich 1995, Hoflich et al. 
1999) either in the photosphere or in the outer scattering 
envelope when the line-of- sight is not coincident with an 
axis of symmetry. Unfortunately, it is not possible from 
polarization measurements alone to decode both the shape 
of the object and the viewing angle and therefore the small 
polarization values reported here could indicate large de- 
parttires from sphericity but viewed close to an axis of 
symmetry, or small departures viewed well away from the 
symmetry axis. Net polarization could also result from 



a spherical envelope in which, as a result of large scale 
clumping, one hemisphere had a surface brightness dif- 
ferent from that of the other hemisphere. Delayed light 
echoes could also be responsible for the observed polariza- 
tion, provided that the time delay is comparable to the 
time scale of the SN luminosity variation. 

We note that no polarization, neither circular nor lin- 
ear, was detected in the radio and this has been inter- 
preted as the signature of a spherically symmetric blast 
wave (Kulkarni et al. 1998). This apparent discrepancy 
might suggest that the radio and the optical radiation were 
generated in regions of different geometry. 

5. HIGH RESOLUTION SPECTROSCOPY 

To search for possible narrow emission lines produced in 

the circumstellar environment of the progenitor and for in- 
terstellar absorption features, we obtained high resolution 
spectra of SN 1998bw around maximum brightness. With 
the exception of a narrow Ha emission line centered at 
2532 km s~^ and probably arising from the underlying HII 
region (sec also Sec. 2.1 and Tinncy et al. 1998), neither 
the usual Nal D absorption lines, which would signal in- 
tervening interstellar dust, nor narrow emission lines were 
detected. However, the signal-to-noise ratio we achieved 
(SNR ~ 20 at Ha and ~ 15 at Nal D) is insufficient to 
exclude the presence of lines, both in emission and ab- 
sorption, with EW< 0.2 A. Using the relation recently 
calibrated by Benetti et al. (in preparation) and assuming 
that the total EW of the Nal D doublet is EW< 0.4 A, one 
can estimate an upper limit for the extinction, Ay <0.2. 
There are two values of Ay available in the literature, i.e. 
Av=Om (Burstein & Heiles 1982) and Ay=0.2 (Schlegel, 
Finkbeiner & Davis 1998). Both estimates are compatible 
with our upper limit, but we will assume Av=0.2 through- 
out this paper due to its more recent determination from 
COBE/DIRBE and IRAS/ISSA maps. 

6. IR SPECTRA: HELIUM AND OTHER ELEMENTS 

The IR spectra are plotted in Fig. 7. The most promi- 
nent feature is the broad emission at 1.08 /xm, which has 
an associated P-Cygni structure with a hint of more than 
one component in absorption. The main component of the 
emission and the absorption minimum at day +8 strongly 
suggest that this feature is due to He I 1.083 /xm with a 
velocity of 18,300+700 km s"-'^. There is nevertheless the 
possibility that transitions from other ions are affecting 
this profile. The absorption does not seem to shift red- 
wards with time as shown in the upper panel of Fig. 8. 
This is to be expected if the helium were confined to an 
outer higher velocity layer either optically thick or thin 
but restricted to a narrow velocity range. A more ex- 
tended layer which had become optically thin on or before 
the date of the first IR spectrum could also in principle 
give rise to this constancy of velocity. The first interpreta- 
tion is reinforced by the fact that at the same phase (i.e. 
22.4 days after the explosion). Si II 6355 A had a much 
smaller velocity (12,000 km s~^) which may more closely 
indicate the photospheric velocity, at least up to day +22 
(see Fig. 6). Since the color temperature of the envelope at 
those phases is not sufficiently high to radiatively or col- 
lisionally excite the required levels of He I, one may con- 
clude that excitation is caused by non-thermal electrons 
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Fig. 7. — The IR spectroscopic evolution of SN 1998bw from +8 to day +51. For presentation the spectra have been vertically shifted by 
an arbitrary amount; +33'* (—0.80, log Fx), +52'' (—1.70). The vertical dotted line is placed at the wavelength of the absorption trough's 
minimum of He I 1.083 //m. The increased noise level in the wavelength ranges 1.35-1.45 /um and 1.8-1.9 /im is due to high atmospheric 
absorption. Spectra are in the host galaxy rest frame. 

Fig. 8. — Upper panel: IR spectroscopic evolution of SN 1998bw from day +8 to day +51 in the region of He I 1.083 /xm. The vertical 
dashed line marks the expansion velocity measured for He I 1.083 /im on day +8. Lower panel: region of He I 2.058 /um at day +51. The 
solid line is a smoothed version of the spectrum. A 350 A wide boxcar filter was used. 

Fic;. 9. — Optical and IR spectra of SN 1998bw at comparable phases. Line identifications from spectral modeling are plotted for the most 
prominent emission features (top) and for the He I lines (bottom). The He marks are placed at the expected absorption positions for an 
expansion velocity of 18,300 km s~^. 



produced by 7-rays coming from the radioactive decay of 
56Co. This was predicted by Chugai (1987) for SN 1987A, 
and shown to be effective by Graham (1988), Lucy (1991) 
and Mazzali & Lucy (1998) in various SNe, after Harkness 
et al. (1987) had identified He I hnes by postulating large 
departures from LTE to explain the abnormal strengths. 

He I has another transition (singlet series) in the near 
IR, namely at 2.058 /zm {2s^S - 2p^P°). Careful inspec- 
tion shows that in all three spectra a weak broad P-Cygni 
feature is present at this wavelength (see Fig. 7). While 
this confirms the presence of He I in the spectra, because 
this spectral region is affected by strong telluric absorp- 
tion, any line strength measurements should be treated 
with some caution. In the lower panel of Fig. 8 the region 
centered on 2.058 /im is plotted for day +51, i.e. when this 
feature appears to be more pronounced. Even though the 
signal-to-noise is poor, a P-Cygni profile is visible and 
from the minimum of the absorption through an expan- 
sion velocity of 13,000 ± 2000 km s~^ is deduced, while 
the measured intensity ratio He 1 1.083/He I 2.058 is about 
30. 

Since these He I lines show P-Cygni profiles associated 
with strong continuum radiation, it is instructive to fol- 
low the qualitative results obtained by an elementary SN 
model in the Sobolev approximation (see for example Jcf- 
fery & Branch 1990). In that framework, the absolute 
emission line intensity depends not only on the line opti- 
cal depth, but it is also proportional to the photospheric 
continuum intensity at the line wavelength. Thus it is not 
surprising that the emission component of He 1 1.083 fxm 
is much greater than that of He I 2.058 /um. There is 
in fact roughly a factor 10 difference in the continuum 
flux at the two wavelengths, and therefore the measured 
intensity ratio translates into an equivalent width ratio 
EW{He 11.083) /EW {He 72.058) -3. The theoretical 
work by Lucy (1991) on the non-thermal excitation of he- 
lium in Type lb SNe predicts a ratio close to one for these 
two He I features, whereas the ratio between the A-values 
is 5.2 (Martin 1987). 

An additional complicating aspect to consider is that the 
He I 1.083 /xm line originates from the metastable triplet 
ground state which can give rise to resonance scattering 
not expected for the He I 2.058 /im line, whose lower level 
is metastable but whose upper level can be depopulated by 
an allowed transition to the singlet ground state. More- 
over, the intensities of the two lines depend also on the 
fraction e of collisional to total depopulations of the upper 
level of the transition. Hence, if one takes also into account 
that the optical depths of the two lines can be different, it 
is not surprising that the observed line ratio deviates from 



the expectations based on Einstein coefficients alone. 

We note that the IR spectra of SN 1987A have shown a 
strong He 1 1.083 /im line starting at day 110 (Meikle et al. 
1989), from which an expansion velocity of 5000 km s^^ 
was deduced. The He I 2.058 iim line was also detected 
with an expansion velocity of 1500 km s~^. Meikle et al. 
(1989) suggested that this is probably due to the fact that 
He I 2.058 /zm was optically thin. For SN 1987A the ratio 
He I 1.083/He I 2.058 in the phase range 110-349 days 
ranged from 23 to 50. In particular, on day 112 it was 
about 35; once one has taken into account the different 
continuum level at the two wavelengths, this turns into an 
equivalent width ratio of about 6. 

Very few other SNe have been observed in the 2 /zm 
region. The Type II SN 1995ad and Type Ic SN 1997B 
have both shown a strong He 1 1.083 /im feature without 
any evidence of the He I 2.058 /iin line (Clocchiatti et al. 
2001). He I 1.083 /zm was identified also in the near IR 
spectra of the Type Iln SN 1998S (Gerardy et al. 2000) 
and in that case, even though very faint, the He I 2.058 /xm 
line was detected at 95 days after V maximum. 

One possible difficulty with the identification of the 
1.083 /im feature as He I comes from the apparent ab- 
sence of the He I lines in the optical spectra at compara- 
ble phases. In Fig. 9 we plot the optical spectra at day 
-|-6, -|-29 and -|-52, i.e. as close in time as possible to the 
IR spectra. At all three epochs, the spectrum is predom- 
inantly an absorption spectrum. Some of the strongest 
features tentatively identified using spectral synthesis are 
marked (top). We also marked the expected position of the 
He I absorption lines if they all had a blue-shift of 18,300 
km s~^ (bottom). The strongest optical line of He I should 
be 5876 A, but this is blended with Na 1 D. He I 5876 A 
is clearly not present as a distinct individual feature at 
18300 km s~^. There is no unambiguous sign of any other 
possible He I optical lines. This does not necessarily con- 
tradict the identification of the 1.083 /xm feature as He I, 
because the optical lines can be suppressed relative to the 
IR ones if non-thermal excitation is at work (Mazzali & 
Lucy 1998). The large expansion velocity of SN 1998bw 
might also dictate greater blending effects in the optical 
region where the number of lines per unit wavelength is 
much larger than in the IR. 

The optical He I lines are at least weaker than in the 
Type Ic SN 19941. This SN also showed the He 1 1.083 /xm 
line (Filippenko et al. 1995) and possibly traces of He I 
contamination in the optical (Clocchiatti et al. 1996). Mil- 
lard et al. (1999) have modeled SN 1994rs spectra and 
could fit the feature near 1.08 /xm using a blend of CI 
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Fig. 10. — Spectroscopic evolution of SN 1998bw from day +29 to day +376. The phases have been computed from B maximum light (May 
10, 1998). For presentation the spectra have been vertically shifted by arbitrary amounts: +45'' (—0.20), +52'' (—0.55), +64'' (—0.95), +73'' 
(-1.20),+94'' (-1.50), +125'' (-1.70), +201'' (-1.90), +337'' (-2.00), +376'' (-2.70). Narrow emission lines visible in the late time spectra 
originate in the parent galaxy. Spectra are in the host galaxy rest frame. 

Fig. 11.— Comparison between SNe 1992A (la, ESO-KP data base), 1998bw, 1996N (lb, SoUerman et al. 1998) and 1996aq (Ic, ESO-KP 
data base, unpublished) at late phases. The vertical dashed line is placed at the rest-frame wavelength of Mg I] 4571 A. 



1.0695 /xm and He I 1.083 /um, the latter being detached 
at 18,000 km s~^. This velocity is similar to what is seen 
in SN 1998bw. However, they were able to obtain an al- 
ternative fit using only Si I lines. Another possibility is 
that the feature at 1.08 /im is due to Mg II 1.091 /um, 
as suggested by Mazzali & Lucy (1998) for the Type la 
SN 1994D. 

Determining the amount of He in the ejeeta is very im- 
portant for constraining the nature of the progenitor and 
its evolutionary stage at the time of the explosion. De- 
tailed spectral modeling, including non thermal effects, is 
necessary for this purpose. Inclusion of helium in the ex- 
ploding progenitor object may even affect our quantitative 
conclusions concerning the ejected mass, mass of radioac- 
tive material, kinetic energy and mixing of the ejecta. 

As the SN ages the 1.083 /im emission feature develops 
a second, redder component (cf. Fig. 7). Possible identifi- 
cations are O I 1.129 /xm and Na I 1.138, 1.140 /xm. The 
latter seems to be possible since its lower level is the up- 
per level of the Nal D line transition which would be well 
populated by resonance scattering. The Na I identification 
implies a velocity close to that of Si II 6355 A at the same 
phase and consistent with the Na I D line absorption. The 
presence of O I 7771 A in the optical is only a weak argu- 
ment in favor of the O I 1.129 /xm identification, since this 
latter line originates from a higher level. 

7. ENTERING THE NEBULAR PHASE 

The spectral evolution of SN 1998bw towards and dur- 
ing the nebular phase is presented in Fig. 10, which shows 
all our spectra taken between day -|-29 and day -1-376. The 
transition from an absorption to an emission spectrum is 
slow and subtle and is marked by a decrease in the con- 
tinuum. 

While the evolution of SN 1998bw in the range 5500- 
9000 A is similar to that of SN Ic events such as SN 1987M 
(Filippenko et al. 1990), the expansion velocities are larger 
and the region between 4000 and 5500 A is dominated (at 
least until about day -1-200) by a wide bump to which 
Fe II transitions probably contribute significantly. The 
spectroscopic differences appear clearly in Fig. 11, where 
SN 1998bw is compared to the Type la SN 1992A, the 
Type lb SN 1996N and the Type Ic SN 1996aq. From 
Fig. 11 the late nebular spectrum of SN 1998bw appears 
to mostly resemble a Type Ic SN, in this case SN 1996aq. 
Mg I] 4571 A and [O I] 6300 A are present in both cases 
but not in the Type la SN 1992A. SNe la in general show 
strong [Fe II] and [Fe III] emission features not appar- 
ent in the same form in the two Type Ic SNe. In fact, 
the strongest spectral resemblance of SN 1998bw at neb- 
ular phases is to the peculiar SN 1985F and the Type Ic 
SN 1987M (Filippenko et al. 1990). Qualitatively the 
appearance of Mg I and O I is consistent with what is ex- 
pected from the nucleosynthesis model of Iwamoto et al. 



(1998), and their absence in models of Type la proposed 
by Nomoto et al. (1997). 

All these facts support the general idea that SN 1998bw 
is related to SNe Ib/c. It might be regarded as an ex- 
treme case among these objects, having large kinetic en- 
ergy, ejecta mass and ejected mass of synthesized ^^Ni, 
while SN 1997ef could represent a less extreme case closer 
in properties to the known SNe Ic (see also Iwamoto et al. 
2000; MazzaU et al. 2000; Stathakis et al. 2000, Branch 
2001). 

At late epochs the most prominent spectral feature is 
the [O I] 6300,6364 A doublet. Fig. 12 (upper and central 
panel) reveals that this line has a complex profile, showing 
an apparently blended emission component blue shifted by 
approximately 2300 km s~^ with respect to the 6300 A 
doublet component. It is present at 201 and 376 days, and 
a similar feature appears to be present in the spectra of 
SN 1996N at comparable phases. This therefore points to 
a common origin in the form of an as yet unidentified line 
(centered at ~6250 A), or a similar structure (a clump?) 
in the line emitting region. 

The [O I] 6300 A line profiles are sharply peaked in the 
nebular phase, as are those of Mg I] 4571 A albeit a lit- 
tle less clearly. Uniform emission throughout the envelope 
should produce a parabolic profile, and therefore we can 
conclude that stronger emission is occurring at lower ve- 
locities in regions closer to the center (see also Maeda et 
al. 2000). A similar effect is apparent in other SNe such 
as the Type Ic SN 1987M and peculiar SN 1985F (Fihp- 
penko et al. 1990) where it is also evident in emission lines 
of other ionic species. 

Since Mg I] 4571 A seems to be the only single line 
in the spectrum with a sufficient signal-to-noise ratio, it 
is of some interest to study its behavior with time (see 
Fig. 13, right panel). At early stages prior to 100 days 
the apparent emission feature must be dominated by some 
as yet unidentified feature because the peak is displaced 
bluewards by 3500 km s~^ and gradually reaches the rest 
wavelength of Mg I] 4571 A at day 125. Prom this phase 
onwards there is an excess in the blue wing causing an 
asymmetry, which becomes more marked as time passes by 
(see also Fig. 10). One explanation for this blue excess af- 
ter day 201 is that it is the remnant of the unknown feature 
that dominated at early phases. It cannot be caused by 
the presence of Mg I] 4562 A which only reaches strengths 
comparable to that of Mg I] 4571 A at densities lower than 
a few 10^ cm^^ and in any case the separation is too small. 
A more interesting and plausible explanation for this blue 
component of Mg I] 4571 A is that it is the counterpart 
of the blue component of [O I] 6300 A mentioned above. 
The velocities and the profiles are similar (Fig. 12, lower 
panel). That being so, it suggests that we are seeing two 
components of the oxygen and magnesium mass distribu- 
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Fig. 12.— Upper and central panel: [O I] 6300,6364 A line profile at two different epochs for SN 1998bw (thick line) and SN 1996N (thin 
line, SoUerman ct al. 1998). The vertical dotted line is placed at the position of the unidentified feature (see text). Velocities are computed 
from the 6300 A doublet component. Lower panel; Comparison between the profiles of [O I] 6300,6364 A (thick line) and MG I] 4571 A (thin 
Une) for SN 1998bw at +376 days. 

Fig. 13.— Evolution of Mg I] 4571 A (right panel) and [O I] 6300,6364 A (left panel) line profiles. 

Fig. 14. — Upper panel: FWHM expansion velocity deduced for Mg I] 4571 A and [O I] 6300,6364 A. Lower panel: evolution of line 
luminosities. Dotted and dashed lines represent linear fittings to the data; line lengths indicate the phase range which has been used for the 
linear regression. The box in the upper right corner of the lower panel shows the region where the flux for [Fe II] 5215 A has been estimated 
and the position of the two continuum levels (dashed and dotted lines) . See text for more details. 



tion in the envelope. 

We estimated the expansion velocity of the Mg nebula 
fitting the Mg I] 4571 A line profile. The results are plot- 
ted in the upper panel of Fig. f4. We emphasize that, 
especially at very late phases, the fitting is not satisfac- 
tory and the values plotted are only indicative. A com- 
parison with other SNe of similar Type is interesting. The 
FWHM of Mg I] 4571 A in the latest spectra of the Type 
lb SN 1983N (225 days, Gaskell et al. 1986) and the Type 
Ic SN 1996aq (270 days) was about 5000±300 km s"^ 
and 7000±800 km s~^ respectively. A similar value (6100 
km s^^) was reported by SoUerman ct al. (1998) for the 
Type lb SN 1996N at about 220 days. These values are 
to be compared with 9800±500 km s~^ measured for the 
same line in SN 1998bw at 201 days past maximum. How- 
ever, this value reduces to 7700±800 km s~^ if we deter- 
mine the velocity using only the red half of the profile. 
Thus it seems that when account is taken of the some- 
what different phases at which measurements were made 
for these SNe, not very significant differences are notice- 
able. 

A similar analysis is more complicated for the [O I] dou- 
blet, both because of its complex profile (cf. Fig. 12) and 
a strong contamination by the Si II 6355 A line on the 
red wing at phases earlier than 3 months (see Fig. 13, left 
panel). We tried to estimate the expansion velocity of 
the O I nebula using the blue wing of the observed profile 
and hence computing the FWHM of the 6300 A compo- 
nent alone. The results are plotted in the upper panel of 
Fig. 14, which shows that the velocities are similar to those 
computed from Mg I]. 

Finally, we have measured the flux of the Mg I] and 
[O I] lines to look for possible deviations from exponen- 
tial decay. The results are plotted in the lower panel of 
Fig. 14, where it appears clearly that at late phases the 
flux in these lines decreases exponentially. The measured 
decline rates are 0.0185 and 0.0170 mag day~^ for the 
two features, respectively. These values are quite simi- 
lar to those measured for the broad band photometry (cf. 
Sec. 8). It should be noted here that Si II 6355 A line 
emission probably dominates the [O I] 6300,6364 A lines 
up to day 94 although there is the possibility that if [O I] 
is optically thick and the expansion velocities significantly 
exceed the doublet separation, most of the [O I] photons 
would be redshifted and confused with those belonging to 
the Si II 6355 A transition. 

The unidentified blue shifted component of Mg I] will 
clearly dominate the photometry of this feature prior to 
day 125 and still contribute afterwards. The fact that 



Mg I] appears to reach the exponential decay line before 
[O I] (cf. Fig. 14, lower panel) could be due to the nature 
of these blending effects. Also other considerations may 
be important such as the fact that the critical density for 
Mg I] is 3 orders of magnitude higher than that for [O I] 
and the expanding envelope results in the former critical 
density being reached well before the latter. Also, since 
Mg I has a much lower ionization energy than O I, emis- 
sion of these different lines may be coming from different 
regions. Finally models of massive stars (Woosley, Langer 
& Weaver 1995) show that magnesium occupies a more 
confined radial distribution than oxygen. This is the case 
not only for extremely massive stars but also for 13-15 Mq 
stars (cf. Fig. 1 of Nomoto et al. 1995). In both circum- 
stances, the region where O exists without Mg is mostly 
He-rich, the He mass fraction being ~0.1. 

Unfortunately, it is not possible to make a similar anal- 
ysis for the presumed Fe lines and particularly the feature 
near 5215 A. In fact even the identification of this feature 
is in doubt. One promising possibility is that it is due to 
the [Fe II] multiplet 19 with additional components due 
to the [Fe II] mufiiplets 16,17,18. [Fe II] mufiiplet 6 may 
then contribute to the blend near Mg I] 4571 A. Never- 
theless, there is no unequivocal evidence for the presence 
of other multiplcts of [Fe II] . Another possibility suggested 
by Filippenko et al. (1990) for SN 1985F and SN 1987M is 
a combination of Fe II multiplets 35,42,49. Contributions 
from [Fe HI] multiplet 1 seems a more remote possibility 
because some important lines are not present. What we 
can state here is that the most prominent feature, the one 
at 5215 A, declines faster than those of Mg I] and [O I]. 
The flux measurement is strongly hampered by the uncer- 
tain position of the continuum and the possible prcscmcc 
of another feature in the red wing (see Fig. 10). To re- 
duce the effect of these problems we computed the fluxes 
by integrating the spectra in the range 5080-5680 A and 
crudely assuming the continuum level. To estimate the 
uncertainty involved we chose the continuum at two ex- 
treme positions. One is the average value in the range 
5680-5700 A, where the fiux drops to a minimum value, 
and the other is the level used in computing the flux of the 
Mg I] 4571 A measurements. 

The results arc shown in the lower panel of Fig. 14, 
where the upper and lower extremes of the error bars rep- 
resent the values obtained using the low and high con- 
tinuum respectively. The conclusion is that the 5215 A 
feature fades at a rate of 0.021±0.001 mag day~^, i.e. 
0.0025±0.001 and 0.004±0.001 mag day-^ faster than [O I] 
and Mg I] respectively. The faster decline of 5215 A is also 
clearly visible, if one compares its intensity with the one 
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of the adjacent Mg I] 4571 A after normalization to its 
peak (sec Fig. 15). This fact might also suggest that the 
5215 A feature does not trace the bulk of the Fe in the 
envelope since one would expect the abundance of Fe to 
be increasing with time as a result of the radioactive de- 
cay of ^^Co. This might be partially compensated by the 
decreasing temperature and density of the ejecta. 

Another point concerns the expansion velocity deduced 
from this feature. Again the measured values are quite 
uncertain, but at day 201 the FWHM velocity is about 
10,900 km s~^ and hence higher than those estimated for 
Mg I] 4571 A (9800 km s-^) and [O I] 6300,6364 A (7600 
km s^^). However, if the 5215 A feature results from a 
blend of multiplct lines of whichever type, this would most 
certainly lead to an overestimate of the expansion velocity 
since the separation of individual lines of the most likely 
multiplets is of the order of thousands of km s~^. Here 
we only note that recent Chandra X ray observations of 
Cassiopeia A, which is consistent with the remnant of a 
massive star explosion (e.g. Fesen et al. 1987), show that 
significant amounts of Fc are present at high velocities, 
comparable to those of lighter elements such as oxygen 
(Hughes et al. 2000). This implies that a relevant fraction 
of the material synthesized in the core has mixed into the 
outer regions during the SN explosion. 

Finally, we have compared SN 1998bw with the Type 
lib SN 1993J and the Type lb SN 1996N at about one 
year after maximum light (Fig. 16). The resemblance of 
the three objects is evident. In particular, SN 1993J and 
SN 1998bw are strikingly similar, the only deviation being 
the broad Ha. which is missing in SN 1998bw; apart from 
this, line ratios and expansion velocities are definitely com- 
parable. SN 1996N also has very similar line widths, while 
the luminosity ratio between [O I] and all other features is 
smaller and the blue [Fe II] bump is less pronounced. 

Despite the great peculiarity shown in the early phases, 
at one year after maximum SN 1998bw is spectroscopically 
indistinguishable from known Type Ib's. Even the high 
expansion velocities measiircd during the first six months 
have slowed down to the values which are typical for other 
Type lb SNe (-5000 km s^^). But the much higher ejected 
mass estimated by the models and the high luminosity, 
which persists also at these advanced phases (cf. next sec- 
tion), confirm that this event was hyper-energetic. 

8. BROAD BAND AND BOLOMETRIC LIGHT CURVES 

The photometric data presented in this paper and those 
published by Galama et al. (1998b), McKenzie & Schaefer 
(1999) and Sollerman et al. (2000) can be collected in a 
unique data set to study the late phase behavior of the 
light curves, as shown in Fig. 17. Fig. 17 shows that the 
SN luminosity follows an exponential luminosity decline up 
to about 300 days after maximum light, while a clear flat- 
tening is present for the latest observed epochs, especially 
in the V band. This is clearly shown in Tab. 4, where we 
have reported the decay rates computed via linear least 
squares fit to the data in the two phase ranges 40-330 
and 300-490 days from B maximum. The values for the 
earliest data are consistent, within the estimated errors, 
with those computed by McKenzie & Schaefer (1999) us- 
ing data in the range 47-171 days {jb = 0.0141 ± 0.0002, 
7y = 0.0184 ± 0.0002, 7/ = 0.0181 ± 0.0002 mag day-^). 



In any case, these slopes deviate clearly from the val- 
ues expected from the ^^Co ^^^Fe radioactive decay in 
the case of complete 7-ray trapping (0.0098 mag day~^), 
as already noticed by McKenzie & Schaefer (1999). On 
the other hand, the measured slopes are similar to those 
reported by Sollerman et al. (1998) for the Type lb 
SN 1996N in the phase range 180 340 days (7^ = 
0.0167 ± 0.0023, 7fl = 0.0172 ± 0.0010, 7/ = 0.0193 ± 
0.0024 mag day~^). These values are even higher than 
those typical for the late phase data of Type la SNe. If 
^^Co decay is powering the late light curve of SN 1998bw, 
there must be a fairly strong leakage in the 7-ray deposi- 
tion, as for SN 1996N (Sollerman et al. 1998). Interest- 
ingly, fast and slow decliners seem to exist within the Type 
Ic class, as pointed out by Clocchiatti & Wheeler (1997). 
The range spanned by late phase decline rates of SNe Ib/c 
is quite broad: jb ranges from 0.006 (SN 1990B, Piemonte 
2001) to 0.022 mag day"! (SN 19901, ESO-KP data base). 
At later phases {t >300 days) a clear flattening of the light 
curve is visible in all passbands, even though the conclu- 
sion about B band is somewhat hampered by the lack of 
data at phases later than 403 days. The decline rate varia- 
tion is 0.0020±0.0013, 0.0079±0.0010, 0.0036±0.0020 and 
0.0043±0.0013 mag day"^ in B, V, R and / respectively. 

In Fig. 18, we present the UVOIR bolometric light curve 
of SN 1998bw, which was constructed using all available 
data by integrating the flux in the optical and near IR 
photometric bands. Magnitudes were converted to fluxes 
using the calibrations of Bessell (1979) and Wilson et al. 
(1972). The measured fluxes were then transformed to 
luminosities by adopting a distance modulus jjL = 32.89 
{d = 37.8 Mpc; = 65 km s'^ Mpc'^), and an ex- 
tinction Kv = 0.2. We emphasize that JHK photometry 
is available only for the early phases and that its contri- 
bution to the total flux is quite significant, ranging from 
~ 25% on day 22.4 (from GRB) to - 35% on day 65.4. 
For later epochs we simply assume that the fractional IR 
contribution remained constant at 35% as for the last mea- 
surement. Clearly this is a major source of uncertainty for 
the late UVOIR bolometric light curve. 

After the maximum peak, where the SN reaches a lu- 
minosity of ~ lO^'' erg s^^, the light curve settles on the 
exponential decay tail starting with day +40, with a de- 
cline rate 7=0.0176±0.0002 mag day-^ Then, after day 
-1-300, the bolometric light curve appears to flatten, grad- 
ually deviating from the slope defined by the data in the 
range 40-200 days (dotted line) so that at -1-490 days the 
difference is about 1.4 mag, which is well outside the er- 
ror bars (cf. Fig. 18). These uncertainties were estimated 
taking into account the presence of possible contamination 
from the unresolved underlying HII region, whose contri- 
bution is difficult to subtract (see below). 

The comparison of the bolometric light curve of 
SN 1998bw with the model of Iwamoto et al. (1998) is 
also presented in Fig. 18. Although the early phases (t< 
60 days) are well reproduced, the observed tail deviates 
from the theoretical prediction. The model must be highly 
energetic to accommodate the early light curve and spec- 
tra. Therefore at later phases 7 -rays should escape quite 
efficiently resulting in a rapid decline. Eventually most 
7-rays escape from the ejecta, and only positrons from the 
^^Co decay contribute to the light curve (the masses of 
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Fig. 15. — The region of Mg I] and [Fe II] at 125, 201 and 376 days (from top to bottom). The spectra have been normalized to the 
Mg I] 4571 A line peak. For presentation, the narrow emission lines emitted by the underlying HII region in the range 4800—5000 A have 
been cut out and the spectra have been smoothed using a 20A wide boxcar filter . 

Fig. 16.— Comparison between Type lib SN 1993J (Barbon et al. 1995, dotted line). Type lb SN 1996N (Sollerman ct al. 1998) and 
SN 1998bw at about one year after maximum light. Spectra have been normalized to the [O I] 6300,6364 A peak and arbitrarily shifted for 
presentation. 



^''Co and ''^Ti in the model are small). Indeed, after day 
~200 the calculated decline becomes slower, and it ap- 
proaches the decay rate of ^^Co around day 400 (Iwamoto 
et al. 1998; Nakamura et al. 1999, 2000a,b). 

Nomoto et al. (2000) and Nakamura et al. (2000b) 
have shown that, while the highly energetic model can ex- 
plain the early light curve, a less energetic model is in 
better agreement with the late light curve. They suggest 
that such a behavior reflects some non-spherical effects in 
the ejecta structure, although the well mixed models of 
Sollerman et al. (2000) do not require asymmetries for a 
reasonable match to the light curves, neither do those of 
Chugai (2000). 

Jeffery's (1999) ad-hoc model, although based on a 
rather different model of the ejecta, also predicts that 
the pure exponential phase should end after about day 
400, and that the light curve should slowly approach the 
^®Co decay slope when eventually only the positrons are 
deposited (provided that other, less abundant radioactive 
species remain unimportant and interaction with circum- 
stellar material does not occur). In this respect it is in- 
teresting to note that a least squares fit to the bolometric 
light curve in the phase range 376-490 days gives a slope 
7 = 0.009 ± 0.001 mag day~^, which would suggest that 
SN 1998bw finally settled on the '"'^Co decay Nevertheless 
a possible alternative is that the flattening is due to the 
unaccounted presence of other sources within the ground- 
based point spread function. In order to reproduce the 
observed behavior, the integrated magnitude of the con- 
taminating objects must be V~22.6. Even though recent 
HST-STIS observations (Fynbo et al. 2000) have shown 
that several objects are present within a radius of 0''5 from 
the SN location (Fig. 1), there is no evidence for such a 
relatively bright source. Hence the conclusion is that the 
observed flattening is, at least to some extent, intrinsic. 
However, this can only be confirmed by follow up obser- 
vations by HST. On the other hand, the effects of possible 
freeze-out, ejecta-wind interaction or even faint echoes 
should not be overlooked. 

As a matter of fact, no signs of ejecta-wind interaction 
are found in the late phase spectra of SN 1998bw (see also 
Sollerman ct al. 2000). Nevertheless, some circumstcllar 
material is expected to be present. All models for SNe Ic 
assume that the progenitor stars undergo intense mass-loss 
and lose their H-He envelopes to become bare CO cores 
before they explode. Thus interaction should be expected 
at some stage. Future observations will be fundamental 
to investigate the circumstcllar environment and the real 
onset of ejecta-wind interaction. 

In Fig. 19 the absolute V light curve of SN 1998bw is 
presented, and compared with the light curves of some 
Type H (1987A, 1979C), lib (1993J), la (1991T, 1992A), 
lb (19901) and Ic (1992ar, 19941) SNe. For SN 1998bw a 
distance modulus /U = 32.89 while the interstellar extinc- 



tion is assumed to be Ay = 0.2 mag (cf. Sec. 5). It is clear 
that the light curve alone would not allow one to classiiy 
this object as a SN Ih/c, nor indeed to distinguish it from 
thermonuclear SNe. 

Even though the internal extinction for SN 19941 is very 
uncertain (here we adopted Ay = 1.2), SN 1998bw was 
much brighter and had a different light curve shape. Un- 
fortunately, late phase data for SN 19941 are not avail- 
able. However, a comparison is possible with the Type 
lb SN 19901. Maximum light was not covered by the ob- 
servations, but the light curve of SN 19901 between day 
+40 and day +120 is very similar to that of SN 1998bw, 
although about 0.5 magnitudes fainter. Aftc;r that the 
luminosity decay rate of SN 19901 was larger than that 
of SN 1998bw. Later on the over luminosity of this SN 
becomes more pronounced. At 350 days the Type lib 
SN 1993J is 1.4 magnitudes fainter then SN 1998bw (see 
Fig. 19). SN 1996N is even less luminous as shown by 
Sollerman et al. (1998), even though the distance modu- 
lus, reddening and epoch of maximum are rather uncer- 
tain. One year after the explosion this SN is about 3 mag- 
nitudes fainter than SN 1998bw. 

As discussed before, we believe that the late time light 
curve of SN 1998bw is powered by the radioactive decay 
of ^'^Co into "''^Fe and therefore the late luminosity can 
be used to constrain the ^^Ni mass ejected by the explo- 
sion. The decay of ^^Co — > ^^Fe releases energy through 
the 7-ray channel in 81% of the cases, while the remaining 
fraction goes into positrons, which annihilate with elec- 
trons producing 7-rays (cf. Colgate & McKee 1969). The 
positron kinetic energy released before the annihilation ac- 
counts for 3.5% of the total ^^Co decay energy (Arnett 
1979, Woosley, Pinto fc Hartmann 1989). 

In the case of SN 1998bw, since the slope of the late 
light curve prior to 300 days is steeper than the expected 
input energy from Co decay, we can assume that in this 
phase range some fraction of the 7-rays escapes thermal- 
ization. At later epochs {t >400 days) the envelope be- 
comes completely transparent to 7-rays whereupon in the 
case of complete deposition of positron kinetic energy (Ax- 
elrod 1980), the late light curve settles on the ^^Co decay 
line. 

Making the crude assumption that the phase of complete 
escape of 7-rays and complete deposition of positrons was 
reached after 400 days and assuming a rate of ^^Co decay 
energy production of s=6. 78x10^ erg g~^ s~^ (Sutherland 
& Wheeler 1984) and a half-life time of 77.12 days (Ar- 
nett 1996), we can estimate the mass of M(^^Ni) using the 
following equation: 

M(56j^i)< ^fj^^ Mo (1) 

^ ' ^ 1.35e-^/i"-26 X 0.035 ° ^ ' 

where L43{t) is the bolometric luminosity in 10*^ erg 
s~^ , f is expressed in days from the explosion and the factor 



Patat et al. 



11 



Fig. 17. — BVRI broad band light curves of SN 1998bw at t >40 days. For presentation the light curves have been shifted by the reported 
amounts. No extinction correction has been applied. The dotted and dashed lines represent least-squares fittings to the data in the ranges 
40-330 and 300-490 days past B maximum. Data are from Galama et al. (1998b), McKenzie & Schaefer (1999), Sollerman et al. (2000) 
and this work. The thick dashed line corresponds to the ®^Co — »^^Fe decay rate, expected for full 7-ray trapping. Epochs refer to the B 
maximum light. 



0.035 is the fraction of total ^^Co decay energy deposited 
by positrons. If we use the four available measurements 
at t > 400 days (Sollerman et al. 2000; this work) and 
average the results, we get M{^^m) < l.Oto^ ^0- The 
errors are by far dominated by the uncertainties in the 
bolometric luminosities and the estimate depends on our 
assumption on the IR contribution at these late phases. 

We note that the model by Nakamura et al. (2000b) 
estimated the mass of ^^Ni to be 0.4 M© from the early 
light curve modeling. This value is not in contradiction 
with what is found here if the SN envelope is not com- 
pletely transparent to 7 rays at t > 400 days. Moreover, 
it must be noted that Nakamura et al. (2000b) have sug- 
gested that positron contribution is not yet dominant in 
powering the light curve at the phases covered by the last 
available observations. 

9. DISCUSSION AND CONCLUSIONS 

As we have shown, SN 1998bw was exceptional in many 
respects, even beyond its possible and probable connec- 
tion with GRB 980425. The luminosity at maximum was 
comparable to that of a SN la but its spectral appearance 
was completely different from that class of objects. On the 
other hand, while the Type Ic classification at maximum 
holds by definition (no H or He, weak Si II), SN 1998bw 
had little in common with objects previously classified as 
Type lb or Ic (cf. also Stathakis et al. 2000), at least 
around maximum. Among all studied Type Ib/c SNe, 
only SN 1992ar might have been brighter than SN 1998bw 
(Clocchiatti et al. 2000a). The only known SN which 
bear some spectroscopic resemblance to this unusual ob- 
ject are SN 1997ef (Garnavich et al. 1997) and SN 1998ey 
(Garnavich, Jha & Kirshner 1998), the former also being 
possibly associated with a GRB (Wang et al. 1998). Even 
so, SN 1997ef was fainter than SN 1998bw, and probably 
produced much less ^^Ni than SN 1998bw (Iwamoto et al. 
1998, 2000). It must be mentioned here that SN 1997cy, 
also conceivably but by no means certainly associated with 
a GRB (Wooslcy ct al. 1999, Germany et al. 2000), was 
extremely bright, probably the brightest SN ever observed, 
but it has shown signs of cjecta-CSM interaction at all 
epochs (Turatto et al. 2000). 

The deviation from any known Type lb SN behavior is 
also noticeable in the early nebular phase, when the spec- 
trum of SN 1998bw is probably dominated by Fe blends in 
the blue and by O and Ca in the red, as already pointed 
out by Patat and Piemonte (1998b). A somewhat similar 
behavior was noticed in SN 1993R by Ruiz-Lapucnte (see 
Filippenko 1997a) and in SN 1990aj by Piemonte (2001). 

The high velocities and the large intrinsic luminosity 
suggest that SN 1998bw was produced by an extremely 
energetic explosion. All published models reach this con- 
clusion, even if the explosion energy, the progenitor mass 
and ejected ^^Ni mass, span quite a wide range (Iwamoto 
et al. 1998, Hofiich et al. 1999, Woosley et al. 1999, 
Nakamura et al. 2000b). This range arises because dif- 



ferent degrees of asymmetry and beaming have been used. 
The exceptionally high value of the kinetic energy of the 
models which give the best fit to both the light curve and 
the spectra of SN 1998bw (3.0 x 10"''^ erg, Iwamoto ct al. 
1998; 5.7 x 10^^ g^g, Nakamura et al. 2000a,b; 6.0 x 10^^ 
erg. Branch 2001) led to the designation of this object as 
a hypernova, even though it does not match the original 
hypernova definition (Paczynski 1998). 

Despite the great individuality shown by SN 1998bw at 
maximum light, one year later its spectrum is very similar 
to that of other Type Ib/c events. This gives qualitative 
support to the idea that the material ejected by SN 1998bw 
was rich in both Fe peak and a— elements. The possibly 
large production of Fe-peak elements makes it important 
to obtain detailed nucleosynthesis calculations and rate es- 
timates of occurrence of such objects, because they might 
have a significant impact on models of galactic chemical 
evolution. Since SN 1998bw was as bright as a SN la, it 
might seem unlikely that such objects would be missed in 
nearby searches. Nevertheless SN 1998bw may have been 
missed had there not been an associated GRB. Does this 
allow us to infer that such kind of explosions must be in- 
trinsically rare at the present epoch? 

But even if they are relatively rare now, their past fre- 
quency may have been considerably higher. SN 1998bw 
is thought to have been generated by a massive star (cf. 
Iwamoto et al. 1998, Woosley et al. 1999) which may 
have been more common at remote epochs, when the rate 
of star formation was higher. Since the time scale for the 
release of iron and oxygen into the ISM by these objects 
is much shorter than for SNe la and significantly less mas- 
sive core-collapse SNe, hypernovae may have played an 
important role in the initial galactic chemical enrichment. 
Future work will concentrate on establishing whether the 
rate of occurrence of such SNe increases with redshift or 
look-back time. It is possible that the discovery of many 
GRB's at significant redshift is already a hint in this direc- 
tion, as is the conclusion for several of them that the light 
curve can be reproduced with a combination of a GRB 
afterglow plus a SN light curve (cf. Bloom et al. 1999). 

This work is based on data collected at ESO-La Silla. 
We express our gratitude to the Visiting Astronomers who 
kindly gave us part of their observing time in order to se- 
cure a good follow-up of this important object. In partic- 
ular, we acknowledge the support we received from Pierre 
Leisy and Alessandro Pizzella during the observations at 
La Silla. We also wish to thank Stephen Holland and Jens 
Hjorth for making the HST STIS image of SN 1998bw 
available to us. Finally, the authors gratefully acknowl- 
edge an anonymous referee for the thorough review, which 
really helped to improve the paper. 
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Fig. 18. — Comparison between the UVOIR bolomctric light curve of SN 1998bw and the hypernova model by Iwamoto et al. (1998, solid 
line). The dotted line represents an extrapolation of a least— squares fit to the data in the range 50-200 days while the dotted-dashed line 
is a fit to the data in the phase range 376—490 days. The dashed line corresponds to the ^®Co — >^^Fe decay rate, expected for full 7-ray 
trapping. The plot in the upper right part of the figure shows the deviations of the observed data from the extrapolation on the early light 
curve (50-200 days, see text). For comparison the bolometric light curve of SN 1987A (Bouchet et al. 1991) is also plotted (filled triangles). 



Fig. 19.— Absolute V fight curve of SN 1998bw compared with SNe 1987A, 1979C (Balinskaya et al. 1980, De Vaucouleurs et al. 1981, 
Barbon et al. 1982), 1993J (lAU Circulars; Lewis et al. 1994; Barbon et al. 1995), 1991T (Phillips et al. 1992, Cappellaxo et al. 1999), 
1992A (Suntzeflf 1996, Cappellaro et al. 1999), 19901 (ESO-KP data base, unpubUshed), 1992ar (Clocchiatti et al. 2000a), 19941 (Tsvetkov & 
Pavlyuk 1995, Richmond et al. 1996). Note that for SN 1992ar the brightest case was chosen (cf. Clocchiatti et al. 2000). The dotted lines 
represent a least-squares fit to the late phase data of SNe 1991T and 1992A. The dashed line corresponds to the ^®Co — >^®Fe decay rate, 
expected for full 7-ray trapping. 



tragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, Cahfornia Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration. 
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oTjanciarci otars 




^^^4UUUUU-t- } 


(days) 


(A) 


r vv xlM 






i.J May yo 


n nn o ^ o 


—9 


oyUu-yzUU 


10 


JN 1 i-rjMMi 


G138-31 


J. 4 May yy 


oU9ob.9 


— 7 


ooou-yuuu 


10 


Danish 


Lii i 737y 


4.d May 9o 


50937.8 


—6 


3400-10250 


20 


Ebu-3.Dm 


/^i— \ J nnoT 

CD32d9yi7 


4.4 May 98 


ouyoY.y 


— 6 


^ A r\r\ Term 


20 


11/00-3.0111-^01. 


rlJJibiUoo 


b.o May yo 


ouyoy.o 


—4 


37oU-7doU 


1 


IN i i-n/L/H 




/.o May yo 


ouy4u.o 


— 3 


Qf;nn nnnn 


10 


Danish 


Lii 1 f o f y 


o.o May yo 


OUy4i.5 


— 2 


QQC^n imr;n 


10 


Danish 


Lii i ( o ( y 


9.3 May 98 


oUy42.o 


— 1 


ooUU-lOioO 


10 


Danish 


Lii 1 737y 


ii.z May 9o 


ouy44.7 


+ 1 


o A nc\ 1 norm 

3450-iU2UU 


10 


Danish 


LI 1 7379 


io.4 May yo 


c;nn Ac n 

ouy4D.y 


+3 


Q A cn 1 nonn 
o4oU-iU^UU 


10 


Danish 


Li i (0(9 


14.4 May 9o 


51)947. 9 


+4 


3350-8900 


10 


Danish 


Lii 7.^79 


16.4 May 98 


50949.9 


+6 


3400-9000 


10 


Danish 


LTT7379 


18.3 May 98 


50951.8 


+8 


9500-25200 


18 


NTT-Sofl 


Hipl06725 


19.4 May 98 


50952.9 


+9 


3250-9000 


10 


Danish 


EG274 


20.3 May 98 


50953.8 


+ 10 


3400-7550 


20 


ESO-3.6m-Pol. 


HD161056 


21.4 May 98 


50954.9 


+ 11 


3400-9000 


10 


Danish 


LTT6248 


22.4 May 98 


50955.9 


+ 12 


3400-9000 


10 


Danish 


LTT7379 


23.3 May 98 


50956.8 


+ 13 


3100-10200 


15 


ESO-1.52 


FeigellO 


29.4 May 98 


50962.9 


+19 


3350-10250 


20 


ESO-3.6m 


LTT7379 


1.3 Jun. 98 


50965.8 


+22 


3400-9000 


10 


Danish 


LTT6248 


8.2 Jun. 98 


50972.7 


+29 


4200-9000 


10 


Danish 


LTT6248 


12.3 Jun. 98 


50976.8 


+33 


9500-25200 


18 


NTT-SofI 


G2720 


24.2 Jun. 98 


50988.7 


+45 


3400-9000 


10 


Danish 


LTT6248 


30.3 Jun. 98 


50994.8 


+51 


9500-25200 


18 


NTT-Sofl 


BS4620,BS6823 


1.2 Jul. 98 


50995.7 


+52 


3400-9000 


10 


Danish 


LTT7987 


13.2 Jul. 98 


51007.7 


+64 


3400-9000 


10 


Danish 


LTT6248 


22.3 Jul. 98 


51016.8 


+73 


3100-10100 


15 


ESO-1.52 


FeigellO 


12.2 Aug. 98 


51037.7 


+94 


3100-10100 


15 


ESO-1.52 


FeigellO 


12.2 Sep. 98 


51068.7 


+ 125 


3400-10200 


20 


ESO-3.6m 


FeigellO 


26.0 Nov. 98 


51144.5 


+201 


3350-10200 


20 


ESO-3.6m 


LTT1020 


12.4 Apr. 99 


51280.9 


+337 


3350-10250 


20 


ESO-3.6m 


LTT6248,LTT7379 


21.2 May 99 


51319.7 


+376 


3350-10250 


17 


ESO-3.6m 


LTT3864, FeigellO 


21.2 May 99 


51319.7 


+376 


4750-6750 


7 


ESO-3.6m 


LTT3864,FeigellO 



Relative to B maximum (JD=2450943.8) 
FWHM of night sky lines 



Note. — NTT-EMMI = ESO-NTT+EMMI+TK2048 (red arm) +TH1024 (blue arm), NTT-ECH = NTT + EMMI echelle mode + TK2048, Danish 
= ESO-Danish 1.54 + DFOSC + Loral/Lesser2048, ESO-3.6m = ESO-3.6m + EFOSC2 + Loral/Lesser2048, ESO-3.6m-Pol. = ESO-3.6m + EFOSC2 
+ Polarimeter + Loral/Lesser2048, NTT-Sofl = ESO-NTT + Sofl + Rockwell Hg:Cd:Tel024, ESO-1.52 = ESC 1.52m + B&C + Ford2048 



Table 2 
ir photometry of sn 1998bw. 



Date 




JD 


Phase'' 


J 


H 


K 


Seeing 


Telescope 






(2400000+) 


(days) 








(arcsec) 




18.3 May 


98 


50951.8 


+8 


13.40±0.04 


13.35±0.04 


13.15±0.04 


0.7 


NTT+SofI 


12.3 Jun. 


98 


50976.8 


+33 


14.41±0.04 


14.25±0.04 


14.20±0.04 


1.2 


NTT+Sofl 


30.3 Jun. 


98 


50994.8 


+51 


15.11±0.04 


14.98±0.04 


14.89±0.04 


0.9 


NTT-l-SofI 



Relative to B maximum (JD=2450943.8). This occurred 14.4 days after GRB 980425. 



Patat et al. 



Table 3 

Late Phase Photometric observations of SN 1998bw. 



Date 


JD 


Phase- 


V 


S 


V 


R 




Seeing 


'i'elescope 




(2400000+) 


(days) 












(arcsec) 




16.1 Mar. 99 


51253.6 


+310 




20.69+0.07 


20.52+0.07 






1.3 


Dutch-0.9 


17.1 Mar. 99 


51254.6 


+311 




20.71 + 0.07 


20.50+0.07 


19.74+0.05 




i.;i 


Dutch-0.9 


8.1 Apr. 99 


51276.6 


+333 




21.10+0.15 


20.69+0.15 


20.09+0.15 




1.1 


ESO-3.6 


12.2 Apr. 99 


51280.7 


+337 


21.13+0.20 








20.03+0.15 


1.1 


ESO-3.6 


21.2 May 99 


51319.7 


+376 




21.59+0.20 


21.43+0.20 


20.83+0.20 


20.61+0.20 


1.2 


ESO-3.6 


17.3 Jun. 99 


51346.8 


+403 


21.69+0.20 


21.91+0.20 


21.70+0.20 


20.87+0.20 


20.76+0.20 


1.0 


ESO-3.6 



Relative to B maximum (JD=2450943.8). This occurred 14.4 days after GRB 980425. 

Note. — Dutch-0.9=0.92m ESO-Dutch + CCD TK512, ESO-3.6=ESO-3.6m + EFOSC2 + Loral/Lesser2048 



Table 4 

Decline rates of SN 1998bw in the phase ranges 40-330 and 300-490 days form B maximum. 



Phase Range** 
N. of points 
7 (mag day-^) 

Phase Range 
N. of points 
7 (mag day-^) 



B 

40.8-311.6 
51 

0.0150±0.0010 

310.6-402.6 

5 

0.0130±0.0008 



40.8-325.6 
57 

0.0180±0.0006 

310.6-489.6 
10 

0.0101±0.0008 



R 

40.8-325.6 
7 

0.0165±0.0019 

311.6-489.6 
9 

0.0129±0.0006 



I 

40.8-325.6 
47 

0.0169±0.0010 

325.6-489.6 

8 

0.0126±0.0008 



Relative to B maximum (JD=2450943.8). 
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